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Abstract 

The 7r + — 7T° mass difference generated by the electroweak interactions of the Standard 
Model is expressed, to lowest order in the chiral expansion and to leading order in a, in 
terms of an integral involving the same correlation function Hlr.(Q 2 ) which governs the 
well known electromagnetic ir + — ir° mass difference. We calculate this contribution 
within the framework of QCD in the limit of a large number of colours N c . We show 
how this calculation, which requires non-trivial contributions from next-to-leading terms 
in the 1/N C expansion, provides an excellent theoretical laboratory for studying issues 
of long- and short-distance matching in calculations of weak matrix elements of four- 
quark operators. The electroweak 7r + — ir° mass difference turns out to be a physical 
observable which is under good analytical control and which should, therefore, be an 
excellent testground of numerical evaluations in lattice-QCD and of model calculations 
in general. 



1 Introduction 



In the chiral limit where the light quark masses are set to zero, the low-energy effective 
Lagrangian of QCD in the presence of electromagnetic interactions has, to lowest order in 
a = e 2 /47r ~ 1/137.03..., an interaction term without derivatives of order O(p ) in the chiral 
expansion jjj 

C eS = --- + e 2 Ctr(Q R UQ L tf) . (1) 



Here, U is the matrix field which collects the octet of pseudoscalar Goldstone fields and 
Qr = Ql = diag[2/3, — 1/3, — 1/3] the right- and left-handed charges associated with the 
electromagnetic couplings of the light quarks. Upon expanding U in powers of the pseu- 
doscalar fields there appear quadratic terms like 

C eS = ■ ■ ■ - 2 e 2 C-L (vr+vr- + K+K~) + • • • , (2) 

Jit 

showing that, in the presence of electromagnetic interactions, the charged pion and kaon fields 
become massive in agreement with the old current algebra result obtained by Dashen Q, 

2e 2 C 

( m K+ ~ m |-o)|EM = (ml+ - m^ )|EM = —P2~ H • (3) 

J It 

In fact, the main contribution to the physical 7r + — ir° mass difference is of electromagnetic 
origin because quark masses do not contribute significantly to the n + — 7r° mass splitting H. 

The effective term in eq. @) results from the integration of a virtual photon in the presence 
of the strong interactions in the chiral limit. It is well known 0-0] that the corresponding 
coupling constant C is given by an integral of a correlation function Hlr{Q 2 ) which is the 
invariant amplitude of the two-point function 

n^(g) = 2i j d A x e **(0|T (^(x)^(O)t) |0> , (4) 

with currents 

L^ = J(x)7^(l- 75 )«(x) and BP = d(x)^(l + j 5 )u(x) . (5) 
In the chiral limit 

KW 2 ) = W - 9^1 2 )^lr(Q 2 ) , (6) 

and 
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C= 8^iJ d Q 2 Q 2l[ LR(Q 2 ), (7) 

with Q 2 the euclidean momentum squared of the virtual photon. The function T1lr(Q 2 ) is, 
for all values of Q 2 , an order parameter of spontaneous chiral symmetry breakdown (S^SB) 
and thus, so is the low-energy constant C. This means that the correlator in eq. (||) has a 
smooth behaviour at short distances, so that the integral in eq. @ converges in the ultraviolet 
region : 

lim Q 4 U LR (Q 2 )^0. (8) 

Q 2 — »oo 
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This behaviour also entails the two Weinberg sum rules || in the chiral limit. It has further- 
more been shown ||[l0| that 

- <2 2 n LR (Q 2 ) > for < Q 2 < oo , (9) 

which in particular ensures the positivity of the integral in eq. ([?]) and thus the stability of the 
QCD vacuum with respect to small perturbations induced by electromagnetic interactions. 

In this letter we shall first show that the same correlation function Hlr(Q 2 ) governs the 
7r + — 7T° mass difference in the presence of the electroweak interactions of the Standard Model. 
This observation and the calculation which follows thereof are of course of phenomenological 
interest per se; but their real relevance lies in the fact that, as we shall show below, they pro- 
vide an excellent theoretical laboratory to study issues of long- and short-distance matching 
in calculations of weak matrix elements of four-quark operators. We also think that this ex- 
ample, which is under good analytical control, provides an excellent testground for numerical 
evaluations in lattice-QCD and for model calculations in general. 



2 The Electroweak tt + — 7r° Mass Difference 

In the presence of the electromagnetic and weak neutral currents of the Electroweak Model, 
the QCD lagrangian £ QC d becomes 

£ Q cd £qcd + qm^r^qR + q^l^L , (10) 

where 

l» = eQ L [A fl (x)-tane w Z fl {x)} + - . - & —T 3 Z^x), (11) 

2 sm 0\y cos 6w 

= eQ R [A^(x) -tan OwZ^x)], (12) 

and 

Qr = Qii = idiag(2,-l,-l), T 3 = diag(l, -1, -1) . (13) 

There are then two possible terms of order 0(p°) in the low-energy effective lagrangian 
generated by the virtual Z integration []: 

+ e 2 tan 2 6 w C z tv(Q R UQ L U^ and - e 2 \ C z tr(Q R UT 3 tf) , (14) 

2 cos z 0\y 

where 

1 3 f°° 2 Q 2 



Cz = ^lJ ^wfMl QY[ ^ Q) - (15) 



These couplings give rise to physical mass terms 

,2 



+ e 2 tan 2 9 w C z %x{QrUQ l U^) = -2e 2 ^^C z ^ (vr+vr- + K+K~) + • • • , (16) 

COS t)\y 



and 



e 



2 1 n^tn _TTrr ttU _ o„2 1 - 1 



C z tT(Q R UT 3 U^ = 2e 2 ^—Cz^(n + TT- + K + K-)+---. (17) 
2 cos z &w cos z &w t4 



lr The charged W-field only couples to left currents and therefore cannot contribute to an O(p ) self-energy. 



2 



Collecting together these weak neutral current contributions with the electromagnetic contri- 
bution in eq. ((2|), we arrive at the quite remarkable result that, to lowest order in e 2 and in 
the chiral limit, 



m 2 K o 



)\sm = (ml+ -mlo)\i 



with 



8^1 ' dQ 



Q 2 



Q 2 + M 2 



2e 2 C s 

fn 



(-Q 2 U LR (Q' 



(18) 



(19) 



In the limit where Mz — ► oo we recover the usual expression of Low et al. [||, while there is 
no induced mass difference in the limit Mz — > where the SU(2) x U(l) gauge symmetry is 
unbroken. This can readily be seen by rewriting the external fields in eq. (|ll|) in terms of the 
usual Wj^ and gauge fields: 



cos 9w 

e 

cos 6w 



Ql - \t. 

QrBJx) . 



B^(x) + -——T 3 W^x) , 
2 sm 0\y 



(20) 
(21) 



Neither the W$ field, which has a pure left-handed coupling to quarks, nor the B^ field, 
which has a left component proportional to the unit matrix, can generate an effective 0(p ) 
interaction of Goldstone fields. When the SU(2) x U(l) gauge symmetry is broken down to 
C/(1) BM , eqs. (18) and ( |l9| ) show that the net effect of the massive physical Z is to reduce just 
a little bit the ir + mass induced by electromagnetism alone. 

It is possible to make a phenomenological evaluation of the integral in eq. ( |T^ ) using the 
fact that the function Hlr{Q 2 ) obeys an unsubtracted dispersion relation 



dt 



1 



t + Q 2 



-ImlWi) - -lmU A (t) 

IT 7T 



(22) 



with ^ImTly(t) and ilmll^i) the physical vector and axial-vector spectral functions which 
are measured in hadronic e + e~ annihilations and in hadronic r decays. The physical spectral 
functions, however, are not quite the "chiral limit" spectral functions which should be inserted 
in eq. (Jl9|) . It is possible to do the appropriate corrections, as discussed e.g. in ref. [11|, but 
for the purposes of this letter we wish to follow a more theoretical procedure. We propose to 



evaluate the integral in eq. (19) within the framework of QCD in the limit of a large number 



of colours N c [12-14], which we shall denote by QCD(oo) for short. 

The spectral function associated with Hlr{Q 2 ) m QCD(oo) consists of the difference of an 
infinite number of narrow vector states and an infinite number of narrow axial-vector states, 
together with the Goldstone pion pole: 



-Imllifl(t) 

7T 



M 2 ) - ]T f A M\5{t - Mi) - f 2 5(t) . 



v 



(23) 



Since Ulr (Q 2 ) obeys the dispersion relation (^ 

-Q 2 U LR (Q 2 ) = f 2 + J2fAM 2 A 

A 



), we find that 

Q 2 



M\ + Q 2 



V 



Q 2 



M 2 +Q 2 



(24) 
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Furthermore, the two Weinberg sum rules that follow from eq. (||) constrain the couplings 
and masses of the narrow states as follows: 

Y,fv M V-Y,f 2 A M l = ti and Efv M V-Ef 2 A M A=0, (25) 
V A V A 

ensuring the convergence of the integral in eq. (|l9|) in QCD(oo) , with the result 



32vr 2 



. A 



Ml 



■775- /aMa log 

m j A s m 2 a 



E 

v 



Ml 



Ml 



Ml 



fvM^log 



Mf 
Ml 



(26) 



The original evaluation by Low et al. Q of the integral in eq. (0) was made by considering 
only the phenomenological contributions from the pion pole and the lowest vector and axial- 
vector states in the narrow width approximation. This approximation, which reproduces very 
well the experimental determination of the physical 7r + — ir° mass difference, can nowadays 
be viewed as the approximation to QCD(oo) which consists in restricting the part of the 
hadronic spectrum which is responsible for SxSB in the channels with the J p quantum 
numbers and 1 + to their lowest energy states. The rest of the infinite number of narrow 
states are then treated as dual to their corresponding QCD(oo) perturbative continuum, and 
therefore do not contribute to order parameters of SxSB. An explicit formulation in terms of 
an effective Lagrangian of this lowest meson dominance (LMD) approximation to QCD(oo) , 
with inclusion of the 0~ and + channels as well, has been recently discussed in ref. |J15| . In 
the LMD approximation to QCD(oo) , the integral in eq. (|i~9| ) gives the result 



a 



SM 



3 fl m 2 



Ml 



m\ Mi _Ma(i 

Ml - Ml g Ml Ml I 8 Ml 



M\ 



M\ 



Ml 



log 



MX 
M 2 



(27) 



where we have expanded in powers of 1 /Ml and retained only the leading term. The first term 
in the r.h.s. of eq. (^7j) is the contribution from the virtual photon integration. The Z-induced 
contribution consists of a "large" log Ml term and a "constant" term. Overall, it represents 
a correction of 0.097% Q to the lowest order electromagnetic contribution. The effect, as 
expected, is very small but nevertheless larger than the present experimental accuracy in 
the determination of the 7r + — tt° mass difference. Although this observation may bear some 
interest by itself, we shall not pursue it any further here. We rather dedicate the rest of this 
note to several theoretically interesting aspects of this calculation. 



3 Ulr(Q 2 ) and the Operator Product Expansion 

As recently discussed in ref. |[(|], besides the two Weinberg sum rules in eqs. (|25|), there are 
further constraints between masses and couplings of resonances on the one hand and the 
local order parameters which govern the operator product expansion (OPE) of the Hlr(Q 2 ) 
function on the other. In particular, it has been shown [TJ] that in the large- N c limit, the 
leading d = 6 order parameter is given by the expression 

lim Q 6 U LR (Q 2 ) = -4 7r 2 (^ + 0(a 2 ))(W) 2 

= J2fvMl-J2flM%, (28) 
V A 

2 For the numerical evaluation, we have taken Mv and Ma as the p(770) and the ai(1260) masses, respec- 
tively. Our normalization convention corresponds to — 92.4 MeV. 
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where (ipip) is the usual single flavour quark bilinear condensate (ip = u, d, or s) in the chiral 
limit, and the Wilson coefficient is the result of a lowest order calculation Jl7| in powers of 
a s . We shall show that it is precisely this d = 6 term which controls the log M§ contribution 
to the Z-induced part of the integral C z in eq. (^). Indeed, following a Wilsonian approach, 
we can split this integral into a low-energy region < Q 2 < fi 2 <C Af| where the Z field 
is integrated out, i.e. where we approximate l/(Q 2 + M^) by 1/M|, and a high-energy 
region /j, 2 < Q 2 < oo where the Z-propagator is fully kept but the Hlr(Q 2 ) function is 



approximated by the leading d = 6 term in eq (^8|). These approximations are equivalent to 
neglecting higher order corrections in 1/M§ and in l//x 2 . We then obtain 



C z 



3 1 
32vr 2 Mf 



32vr 2 Mf 



E/iMilog^ 



E/^iog^ 



V 



V 



log^, 



(29) 



where the first line is the result from the low-energy region and the second line the one from 
the high-energy region. We observe that the separation scale fi 2 cancels in the sum of the two 
terms; in other words, there is an exact matching between the long-distance contribution and 
the short-distance contribution. This cancellation also occurs in the LMD approximation to 
QCD(oo) and the result coincides then with the one given in eq. (^). The reason why we can 
exhibit this exact matching is of course due to the fact that the function Hlr,(Q 2 ) is an order 
parameter for all values of Q 2 and therefore has duality properties under the transformation 
Q 2 1/Q 2 which, in the large- iV c limit, we have been able to work out explicitly. 

It is also interesting to look at eq. (|29| ) from the point of view of effective field theory. 
The full integral in eq. ( |l5| ) can be split as follows 



C z = lim {C Z (M Z -> oo, f i) + SCz(M z ,^)} 



/i— >oo 



(30) 



where 5Cz{Mz, /u) = Cz — Cz(M z — » oo,/x). Each of the two terms in the sum is now 
UV-divergent and this is why an ultraviolet cutoff \x is needed. The first term in eq. ( |30|) 
corresponds to the result from the effective low-energy theory, where the Z field has been 



integrated out and it is precisely given by the first line in the r.h.s. of eq. (29). The second 
term 5C Z (M Z , fi) corresponds to the matching condition in the effective field theory language 
and the result is precisely the one given in the second line of eq. (|29|). Formally one can 
set \x = M z (the effective field theory scale in this case) in these expressions since they are 
actually n independent. Then 5C Z (M Z ,M Z ) = and C z = C Z (M Z -> oo,M z ). 

We are now in a position to discuss how the above calculation proceeds in the framework 
which is usually adopted to tackle weak matrix element calculations, fi 



4 Four— Quark Operators 

The relevant term in the Lagrangian of the Standard Model which is responsible for the 
Z-induced contribution to the tt + — ir° mass difference is the neutral current interaction term 



C 



KC 



2 sin 9\v cos 6 



q L ^T z q L - 2 sin 2 9 w q L ^Q L q L - 2 sin 2 6 w q R ^Q R q R 



(31) 



See e.g. the lectures of A. Buras in ref. [IL8| and references therein. 
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When looking for the induced effective Lagrangian of order O(p ) which contributes to Gold- 
stone boson masses, it is sufficient to consider left-right operators. In the absence of the 
strong interactions, the effective four-quark Hamiltonian which emerges after integrating out 
the Z field is 



-1 e 2 
Ml cos 2 9 W 



sin 2 O w Qlr ~ ^ (qL7f,T 3 q L ) {qR^QRqR) 



where 



Qlr = {qLi^Qiqi) (qR^QRqR) , (33) 

and summation over quark colour indices within brackets is understood. In fact, to O(p ), 
only the first term proportional to the four-quark operator Qlr can contribute. In the 
presence of the strong interactions, the evolution of Qlr from the scale M§ down to a scale 
/j? can be calculated in the usual way, provided this fj? is still large enough for a perturbative 
QCD (pQCD) evaluation to be valid. In the leading logarithmic approximation in pQCD, 
and to leading non-trivial order in the 1/N C expansion, the relevant mixing in this evolution 
which we need to retain is simply given by 

Qlr{M\) = Q L r^ 2 ) log ^fD RL (^) + • • • , (34) 

where Drl denotes the four-quark density-density operator 

Drl = e q e q'(l' 'lQr)(QRQl) > ( 35 ) 
<?,<?' 

with e q and e q > the quark charges in units of the electric charge. This can be seen as follows: 
in the MS renormalization scheme, the full evolution of the Wilson coefficients cq and cd of 
the operators Qlr and Dlr at the one loop level is governed by the equations (subleading 
contributions in the 1/N C expansion have been neglected) 





-3 + - 




(36) 



with boundary conditions: cq(Mz) = 1 and cd(Mz) = 5^. The result in eq. Q34| ) follows 
when taking cd(Mz) = 0, which is appropriate when keeping the one-loop leading log only, 
and from the off-diagonal term in the (transposed) anomalous dimension matrix. 

We are then confronted with a typical problem of bosonization of four-quark operators. 
The bosonization of Drl is only needed to leading order in the 1/N C expansion. To that order 
and to order O(p ) in the chiral expansion it can be readily obtained from the bosonization 
of the factorized density currents, with the result 

D RL = £ e q e ql (q' L q R )(q R q' L ) -> 2B^f x 2£^tr (UQ L U^Q R ) , (37) 
g,q' 



4 See e.g. the lectures in ref. Jt^] and references therein. 
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where B is the low energy constant which describes the bilinear quark condensate in the chiral 
limit, 

ft 

We find that the overall contribution of the term proportional to the Drl([jl 2 ) four-quark 
operator, which we denote Cz\d rl , is given by the expression 



B = . (38) 



c ^»=4 3 ?i w>2 ^ los f ■ (39) 



and it is exactly the same result as the one coming from the d = 6 term of the OPE in the 
previous calculation, the second line in eq. (^), once the duality constraint in eq. ( f28| ) is taken 
into account, i.e. precisely the same contribution as required by the matching condition in the 



effective field theory analysis above. Equation ( p9| ) does not change in an MS renormalization 
scheme when restricted to the one-loop leading log. 

The problem is then reduced to the bosonization of the operator Qlr{^ 2 )- We are con- 
fronted here with a typical calculation of a hadronic matrix element of a four-quark operator, 
in this case the matrix element {ir + \Q lr{^ 2 )\^ + ) ■ The factorized component of the opera- 
tor Qlr, which is leading in 1/N C , cannot contribute to the 0(p ) term of the low-energy 
effective lagrangian. The contribution we want from this matrix element is therefore the 
next -to-leading one in the 1/N C expansion and it requires the evaluation of non-factorizable 
four-quark matrix elements which is a priori a highly non-trivial task. Yet, we know that 
there is a straightforward integral representation of the observable we want to compute in 
terms of a two-point function, and we have succeeded in doing the calculation that way 
within the l/iV c — expansion framework. This seems to indicate that there should be a way 
to do the calculation of the matrix element {^ + \Qlr{^ 2 )\^ + ) as well. Indeed, we shall next 
show that, for this particular matrix element, our present knowledge of analytic results in 
non-perturbative QCD allows us to do the calculation exactly to next-to-leading order in 
the 1/N C expansion. The deep reason behind it is that we know from the discussions in the 
previous sections that only two-point functions can appear in the dynamics of this matrix el- 
ement, and therefore it is enough to have a formulation of the low-energy effective lagrangian 
of QCD(oo) compatible with the OPE constraints at short distances for two-point functions. 
These are precisely the constraints which have been recently discussed in ref. [Oj. 

The calculation proceeds along much the same lines as first suggested in papers by 
Bardeen, Buras and Gerard [2^,21] sometime ago 5 , except that we shall go beyond loops 



generated by Goldstone particle interactions alone in order to achieve a correct matching with 
the logarithmic scale dependence of the short-distance contribution in eq. (j39|). We have eval- 
uated the matrix element (7r + 1 (/x 2 ) 1 7r + ) within the framework of an effective Lagrangian 
which is a straightforward generalization to an arbitrary number of massive J p = 1~ and 
J p = 1 + mesonic states of the effective Lagrangian corresponding to the LMD approximation 



to QCD(oo) recently discussed in ref. [15|. Keeping only the terms that are relevant for the 
evaluation of the 0(p ) contributions, we have 

C eS = JftrD^U + D»U + L w trU + F£ u UF LjlMU 



See also refs. [E2H241 and references therein for more recent work. 
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V 

E /vtrV/f - ^ E /AtrV/^ + - 



(40) 



with the same notation as in ref. [15], Lio being one of the 0(p ) constants of the Gasser- 
Leutwyler Lagrangian [25]. 

The bosonized expressions of the currents (qLl^Qiqi) and (qR'j^QRqR) are obtained from 
the resulting effective action T e fj at tree level, 



(qL-fQLqL) - (Ql) 



sr. 



cff 



I.I 



61% (x) 



and (qR^QaqR) -> {Qn)ij 



ST 



cff 



<5r^ J (a;) 



(41) 



where i, j = u,d,s are flavour indices. The contributions from a loop of virtual pions, vector 
and axial-vector resonances are computed by applying a finite ultraviolet cut-off fi 2 to the 
integral over the corresponding (euclidean) virtual loop momentum. One technical point, 
as already discussed in ref p7| (see eqs. (7.24) to (7.26) of this reference), requires however 
attention. It is the fact that to next-to-leading order in the 1/N C expansion, and besides the 
loops generated by the bosonized factorized currents that we have just discussed, there are 
also contact terms generated by the second variation of the effective action. In our case these 
terms are: 



^contact 



s 2 r. 



iQLUQR ^( x )Sr^ {x) 
{-2fl tr (QrUQlU*) 6(0) - QL W tr (QrUQ l W) DS(0) 



+ 



•} , (42) 



where the J-function contributions are to be interpreted within the same cut-off regulariza- 
tion : 



S(0) 



7T 



" dQ 2 Q 2 = 



and 



□5(0) 



7T 







(43) 



Adding the two types of contributions, the final result reads 

+1 



(n+lQLRifi 2 )^ 



16vr 2 
3 



3 4 

o/ 2 



+ 



f" 2 

J 7T 



dQ 2 Q t 



E 

v 



fv 



Q 2 + M 2 



E 



Q 2 + M 2 



(44) 



The contribution of the Goldstone bosons alone corresponds to the two terms in the first 
line of the r.h.s. of eq. (44). They display a typical polynomial dependence with respect 
to the cut-off fi, which can hardly provide a reasonably good matching with the logarithmic 
scale dependence coming from the short-distance contributions in Cz\d rl - In fact, in an 
MS regularization scheme, as commonly chosen for the evaluation of the short-distance 
Wilson coefficients, these power divergences will automatically disappear. Simply adding 
higher resonances does not by itself solve the problem of matching the long- and the short- 
distances either; however, when the information coming from the short-distance properties 
of the correlator in eq. 



), i.e. the two Weinberg sum rules fl25[) and the sum rule 
-4Lio = X;/v-E/l> 

V A 



(45) 
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is taken into account, the result of eq. (H) can indeed be recast into a form which reproduces 
the first line of eq. (29). Notice that in an MS regularization scheme, the integral in eq. 
ij) should be understood in n = 4 — e dimensions and therefore multiplied by ji A ~ n with \i 



the MS regularization scale. The result, when combined with eq. fl39j ) finally yields eq. ( |29[) 
once again. We insist on the fact that, regardless of the regularization one chooses, the 
calculation we have done of the matrix element (tt + \Qiji(ij 1 2 )\tt + ) is an exact calculation to 
next -to-leading order in the 1/N C expansion and to O(p ) in the chiral expansion. 



5 Comments and Outlook 

The main purpose of the above analysis has been to show that the calculation of the elec- 
troweak contribution to the 7r + — 7r° mass difference in the chiral limit is an interesting theo- 
retical laboratory for testing issues connected with the matching of long and short distances 
in the evaluation of hadronic weak matrix elements in general. 

On the one hand, upon identifying the relevant low-energy constant in terms of the 
appropriate QCD correlator, we have been able to obtain, in the large-iV c limit, an exact 
result in terms of resonance parameters with masses and couplings constrained by the short- 
distance properties of this correlator. 

On the other hand, due to the presence of the large scale set by Mz, we were able to 
proceed as done usually in the study of weak non-leptonic processes. In this respect, we en- 
counter a typical situation: the effective Hamiltonian at tree level is given by a current -current 
four quark operator Qlr, which, upon taking into account short-distance pQCD corrections, 
mixes with the density-density operator Drl- The matrix element of the latter to O(p ), 
which is only needed at the leading order in the large- iV c limit, is easily computed. In the 
same limit, (7r + |Q_L_R(/i 2 )|7r + ) vanishes, and next-to-leading, non-factorizable, contributions 
have to be included. A suitable matching of the long-distance contribution with the loga- 
rithmic scale dependence coming from the short-distance contribution cannot be achieved, in 
this case, by evaluating the relevant hadronic matrix elements in terms of an effective theory 
incorporating only the Goldstone bosons. When the contribution of vector and axial-vector 
states satisfying the appropriate short-distance constraints are included, we find an exact 
matching. 

We think that the observations which follow from this example point towards promising 
perspectives for a systematic determination of hadronic weak matrix elements within the 
framework of QCD(oo) . The required steps are the following: first one needs an identification 
of the low-energy constants of the AS = 1 and AS = 2 effective chiral Lagrangian in terms of 
integrals of QCD Green's functions, i.e. the equivalent of our eq. (f[9|) above. Next one should 
proceed to the study of the constraints which the OPE imposes on these Green's functions 
in the large- iV c limit, i.e. the equivalent of the Weinberg sum rules in eqs. ( p5| ) and the sum 
rules in eqs. ( |28|) and (|45|). The final step is the construction of an effective Lagrangian, along 
the lines shown in ref. flj~5| , which incorporates these constraints. This program is at present 
under study. 

In the meantime, it would be interesting to see a calculation done using numerical simula- 
tions of lattice QCD of the matrix element (vr + |Qx/j(// 2 )|7r + ), so as to be able to appreciate the 
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performance of the lattice approach in a simple case which we understand well analytically. 
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